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We study non-Markovianity as backflow of information in two-qubit systems. We consider a 
setting where, by changing the distance between the qubits, one can interpolate between independent 
reservoir and common reservoir scenarios. We demonstrate that non-Markovianity can be induced 
by the common reservoir and single out the physical origin of this phenomenon. We show that two- 
qubit non-Markovianity coincides with instances of non-divisibility of the corresponding dynamical 
map, and discuss the pair of states maximising information flowback. We also discuss the issue of 
additivity for the measure we use and in doing so, give an indication of its usefulness as a resource 
for multipartite quantum systems. 
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I. INTRODUCTION 

Exploiting quantum correlations, such as entangle- 
ment, as a resource for quantum information revolves 
around the major obstacle of preserving them from the 
detrimental interaction with the environment [TJ [5] . Ef- 
ficient schemes for the protection of correlations gener- 
ally require detailed knowledge of the microscopic pro- 
cesses that induce environmental noise. To this aim we 
consider a simple bipartite system of two spatially sep- 
arated qubits immersed in a non-Markovian structured 
environment. In such an environment memory effects ex- 
ist as a result of complex system-environment coupling, 
and consequently quantum properties can be temporar- 
ily restored. Motivated by the future prospect of scal- 
able quantum devices requiring multipartite systems, and 
in view of the fact that non-Markovianity is a resource 
for quantum technologies [3] , we study how environment- 
mediated interactions enhance memory-keeping proper- 
ties and affect the behaviour of non-Markovianity in two- 
qubit channels. Our main focus is on correlated noise 
induced by the common environment. By changing the 
separation between the qubits we are able to interpolate 
between local and common reservoir scenarios and ex- 
plicitly study the effect of environment mediated correla- 
tions on the degree of non-Markovianity of the two-qubit 
dynamical process. 

In order to systematically investigate the ability of 
an open system to regain quantum properties, it is 
important to unambiguously define and quantify non- 
Markovian dynamical behavior. Here we consider two 
measures of non-Markovianity recently introduced by 
Breuer, Laine and Piilo (BLP) [3] and Rivas, Huelga and 
Plenio (RHP) [5j. These measures are based on simi- 
lar but subtly different definitions of non-Markovianity, 
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i.e., backflow of information from the environment to the 
system and non-divisibility of the dynamical process, re- 
spectively. The definitions agree for some dynamical pro- 
cesses, however, this is generally not the case [SHE]- The 
origin and extent of differences of these measures is an 
active area of research, albeit with most current stud- 
ies focusing on single qubit dynamical processes. The 
present study, comparing and contrasting the measures 
for a two-qubit system, is therefore also of fundamental 
interest. 

Non-Markovianity measures characterize and quan- 
tify certain properties of quantum dynamical maps and, 
therefore, do not depend on the initial state of the open 
quantum system. In practice this means that they are 
defined via optimization over all possible initial states. 
As a consequence they are increasingly difficult to calcu- 
late for multi-qubit systems. Moreover, these measures 
are in general not additive, hence it is not possible to 
reduce the calculation of non-Markovianity of an n-qubit 
channel to the one of a single qubit channel even when 
the qubits are subjected to identical local uncorrelated 
reservoirs [5]. All previous studies on non-Markovianity 
as information flow in two-qubit systems focus on the 
simpler case of independent environments [5Vll2j. The 
model considered here is, to the best of our knowledge, 
the first study of BLP and RHP non-Markovianity in a 
common environment scenario, demonstrating how en- 
vironment mediated correlations create non-Markovian 
dynamics via correlated noise. 

This article is organized as follows. In Sec. II we in- 
troduce our model and the master equation describing 
the dynamical evolution of the system. In Sec. Ill and 
Sec. IV we study the two non-Markovianity measures 
mentioned above in relation to relevant system parame- 
ters, commenting on the optimizing pair and the additiv- 
ity property of the BLP measure, and finally and in Sec. 
V we summarize our results and present our conclusions. 
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II. THE MODEL 

We consider a system of two qubits, spatially separated 
by distance 2D and dephasing under the influence of a 
bosonic reservoir with a non-trivial spectral density func- 
tion. The model has been introduced originally in Ref. 
[13 HI] j however we will later specify our calculations on 
a physically realistic system with a more complex spec- 
tral density function [T5]. When D — > oo the dynamical 
map factorizes into the product of the dynamical maps 
of the individual qubits, = $A(i) <8> ®B(t), where 

subscript A refers to the first qubit and B to the second 
one. For finite distances D, however, this is not the case 
and the qubits are subjected to correlated noise, i.e., to 
environment mediated interactions. On the one hand, 
the conditions for this system to induce non-Markovian 
dynamics on a single qubit have been recently studied in 
Refs. [ini E] ■ On the other hand, a link between the 



dynamical evolution of the qubits and their spatial sepa- 
ration has been considered in Refs. [HI [19] , even if not in 
connection with reservoir memory effects. Our goal is to 
study the interplay between the environmental spectrum 
and the spatial separation in creating non-Markovian dy- 
namics for this two-qubit model, focusing especially on 
the ability of correlated quantum noise to induce a back- 
flow of information. 

We have derived the master equation describing the 
dynamics of the two qubits starting from the exact an- 
alytical expression of the density matrix of two qubits 
at all times given in Refs. [131 [T3] . For the model here 
studied it is possible to solve the dynamics of the to- 
tal closed system exactly and obtain the reduced system 
dynamics by taking the partial trace. Following straight- 
forward calculations, the master equation can be recast 
in the following Lindblad-type form with time dependent 
coefficients. 



dp = li{t) - 72 0) 
dt 2 

7i (*)+ 72(f) 



0-J{(^-"?)(o?-of),p} 



(1) 



where the decay rates are given as the sum and the dif- 
ference of 



71 W = 



l rfTo(t) 
2 dt ' 



72 (t) = 



1 dS(t) 
4 dt ' 



(2) 



Here To (t) is the decoherence function describing individ- 
ual dephasing of each qubit and S(t) encapsulates the cor- 
rections to the single qubit dephasing arising from envi- 
ronment mediated interactions. In the limit of the qubit 
separation going to infinity the latter vanishes and the 
master equation reduces to a sum of two Lindblad-type 
terms describing the individual dephasing of each qubit 
with decay rate 71(f). We note that the single qubit de- 
coherence function To(t) and the correction term S(t) are 
connected to the collective decoherence functions, 



T±(t)=2T (t)±5(t), 



(3) 



of Ref. [T3j . Physically T + (t) and r_(t) describe the 
decay of the super- and sub-decoherent two-qubit states, 
|*±) - ^(|10) ± |01)) and |*±) = i(|00) ± |11», re- 
spectively. Decay of the super-decoherent state is en- 
hanced by environment mediated interaction, whereas 
the decay of the sub-decoherent state is suppressed by 
it. It should be noted that the explicit form and physical 
meaning of the decays rates r±(£) appearing in the two 
qubit channel remain general while their functional form 
is dependent on the specific model studied. 
Initially the two qubits dephase individually as if cou- 
pled to independent environments. At some later time, 



if the qubits are sufficiently close together, they are able 
to "see" each other via correlations arising through an 
interaction mediated between them by the environment, 
formally contained in 5(t). This interaction becomes in- 
creasingly weak as the separation between the qubits in- 
creases, effectively approaching zero at sufficiently large 
distances [T5]. This model is therefore particularly in- 
teresting because it allows insight into these two possi- 
ble physical scenarios through simply changing the sep- 
aration between the qubits. In the following Sections we 
study the dynamics of the decay rates, connecting them 
to two prevailing definitions of non-Markovianity. 



III. DECAY RATES AND DIVISIBILITY 

Useful information on the dynamical properties of the 
two-qubit system is contained in the specific time de- 
pendence of the decoherence rates appearing in the mas- 
ter equation 0. A completely positive, trace preserv- 
ing (CPT) map $(t,0) : p(0) n- p(t) is called divisi- 
ble if $(t,0) = $(t,t)®(t',0), where ®{t,t') and $(t',0) 
are CPT maps Vf'. For master equations in the diag- 
onal form, such as the one of Eq. |l]), the divisibil- 
ity property is reflected in the coefficients of the mas- 
ter equation. Divisibility, in turn, is crucially related to 
non-Markovianity. For all known non-Markovianity mea- 
sures Af a finite value Af ^ implies that the dynamical 
map is non-divisible. In particular, the RHP measure ex- 
plicitly quantifies non-divisibility of the dynamical map, 
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Nrhp if and only if the dynamical map is non- 
divisible [5]. 

For time-local master equations in Lindblad-type form 
the map is divisible whenever all time-dependent deco- 
herence rates, i.e., 71 (t) ± 72 (t) in the case of our model, 
are positive at all times. Conversely, if at least one of 
the time-dependent coefficients temporarily takes nega- 
tive values, then the dynamical map is not divisible. 

Since the negativity of the decay rates depends cru- 
cially on the spectral density function characterizing the 
properties of the environment and the qubit-environmcnt 
coupling, in the following, for the sake of concreteness, we 
will focus on a specific experimentally realizable physical 
system. The set-up we consider comprises of two im- 
purity atoms dephasing due to an interaction with an 
ultracold bosonic Rubidium gas in a Bose-Einstein con- 
densed (BEC) state, introduced in detail in Appendix A. 
Each atom is trapped in double well potential, forming a 
qubit system where the two qubit states are represented 
by the occupation of the atom in either the left or the 
right well. The two impurity atoms are separated by dis- 
tance 2D > 8L, with 2L the distance between the two 
minima for the double-well potential, i.e. the "size" of the 
qubit. It has been shown in Refs. [TBI [T7] that this qubit 
system is particularly sensitive to non-Markovian effects. 
The properties of the bosonic gas can be controlled by 
modifying its trapping potential (in this way one may 
change the effective dimensionality of the BEC) and/or 
changing the boson-boson scattering length as from it's 
natural value a^b via Feshbach resonances (allowing in- 
terpolation from a free to an interacting background gas). 
The dynamics of a single impurity immersed in the BEC 
gas was studied in Ref. [16] . where it was reported that 
the dynamics of the single dephasing qubit is sensitive to 
the scattering length as of the BEC. A critical value of 
the scattering length, denoted et^ 4t and depending on the 
dimensionality of the BEC, was found to specify where 
the crossover between Markovian and non-Markovian be- 
havior occurred. If the background gas is non-interacting 
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or very weakly interacting, i.e. < as < ofj 



the 



qubit dynamics is Markovian. If the background gas has 
stronger interactions, i.e. aB>af qubit dynamics has 
signatures of reservoir memory effects. 

We now study the non-divisibility of the two-qubit dy- 
namical map, signalled by temporarily negative decay 
rates. Choosing the BEC scattering length and the dis- 
tance between the two qubits appropriately, we can con- 
sider the cases when the single qubit dynamics is Marko- 
vian or non-Markovian and the cases when the qubits are 
effectively in local or common environments. 

We begin by looking at the case in which D/L ^> 1, 
i.e., the qubits are subjected to uncorrelated noise, and 
values of the scattering length such that each one-qubit 
dynamical map is Markovian, i.e., 71 (t) is always positive 
(See Fig. 1, upper row). We note that, while 72W ~ 
at all times as expected, 71 (i) increases initially but, 
after a certain time, decreases to zero meaning that the 
dephasing of the individual qubits stops. 
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FIG. 1. (Color online) Time-dependent coefficients 71 (t) (blue 
dashed line) and 72 (£) (blue solid line) are shown in the first 
column, with 71 (i) + 72 (i) (red dashed line) and 71 (t) — 72 (t) 
(red solid line) in the second column. The scattering length is 
chosen to be 0.02anb < a c g l1 to ensure the one qubit channel 
is Markovian for all times t and am, is the natural scattering 
length of a Rubidium condensate, see Appendix A. The insets 
zoom in on small values of the decay rates. 



Keeping all other parameters fixed, we now decrease 
the distance between the qubits and we see that, after a 
short time interval, the cross-talk term 72 (t) begins tak- 
ing negative values (See Fig. 1, lower row). More pre- 
cisely, the maximum negative value of 72(0 is reached 
when 71 (t) is already significantly small. As a conse- 
quence the decoherence rate 71 (i) + 72 (t) takes negative 
values indicating that the two-qubit dynamical map is 
non-divisible. This feature is not specific to the values 
of parameters chosen in Fig. 1, but occurs for any ini- 
tial values of such that, when D/L 3> 1, the map is 
Markovian. Therefore we can uniquely link the presence 
of correlated noise to the violation of divisibility and, 
hence, to a non-zero value of the RHP non-Markovianity 
measure. 

To conclude, the presence of the "cross-talk" term 72 (i) 
is sufficient to break divisibility in the combined system 
regardless of the fact that the individual qubit dynam- 
ics is divisible. More generally, whenever the qubits are 
subjected to Markovian local dynamics, i.e. 71 (t) (a term 
independent of qubit separation) is always positive, non- 
divisible non-local dynamics of the two-qubits can result 
only from the action of correlated noise. 

In order to gain further insight in the microscopic phys- 
ical processes of the composite system, we will look in the 
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next section to the behavior of information flux. Our aim 
is to see whether the cross-talk term acts as the reservoir 
memory, i.e., if it is responsible for backflow of informa- 
tion. To answer to this question we need to study the 
BLP non-Markovianity measure which, as we mentioned 
in the Introduction, does not always coincide with the 
RHP measure. 



IV. INFORMATION FLUX AND 
NON-MARKOVIANITY MEASURE 

The BLP measure of non-Markovianity is based on the 
behavior of information flux between system and envi- 
ronment [4]. According to this definition, a system is 
non-Markovian when part of the information that was 
lost into the environment due to decoherence and/or dis- 
sipation is temporarily restored in the system. To quan- 
tify this process one can use the trace distance between 
two states as a measure of their distinguishability and 
equate decrease of distinguishability with loss of infor- 
mation about the system. Hence, the rate of change in 
trace distance D(pi(t), p2(t)) = \ Tr | p\ (t) — p2(t)\, where 
Pi,2{t) = ${t, 0)pi ) 2(0) are two initial states evolving un- 
der the dynamical map <£(£, 0), can be interpreted as an 
information flux a, 

a(t, Pl ,2(0)) = j t D( Pl (t),p 2 (t)). (4) 

A decrease in the trace distance represents information 
flowing from the system into the environment and an 
increase represents a back flow of information from the 
environment to the system. The total increase of dis- 
tinguishability over the whole time evolution quantifies 
the total amount of information flowing from the envi- 
ronment back into the system so therefore is defined as 
a measure for non-Markovianity: 

Nblp = max / dsa(s). (5) 

The time integration is extended over all time inter- 
vals for which a is positive and the maximum is taken 
over all pairs of initial states, pi,2- Therefore the non- 
Markovianity is defined as the maximal amount of in- 
formation that the system can possibly recover from its 
environment. A dynamical map which is divisible always 
describes Markovian dynamics in terms of information 
flow, however the converse is generally not always true 

[SHE]. 

The maximization in Eq. ^ complicates the calcula- 
tion of the BLP measure of non-Markovianity. In Ref. 
|20j the mathematical and physical properties of the op- 
timal pairs corresponding to the maximal backflow of in- 
formation from the environment to the system are char- 
acterized, simplifying the maximization, but finding the 
optimal pair for a two-qubit dynamical process still re- 
mains a difficult task. To this aim, we calculate the mea- 



sure by randomly generating pairs of states and classi- 
fying them into pure, mixed, separable and maximally 
entangled states including the Bell states. 



A. Non-Markovianity Dependence on Qubit 
Separation 

We will now study how Nblp changes with the dis- 
tance between the qubits. We first consider the case of a 
weakly interacting BEC reservoir inducing local Marko- 
vian dephasing of the single qubits (Fig. [2j left plot). 
The plot clearly shows that the non-Markovianity mea- 
sure decreases as the qubits become further and further 
apart and vanishes for a long enough spatial separation 
between the qubits. This is consistent with the previous 
results (see Fig. [I] first row) : when the qubits are very 
far apart, D/L > 30, they effectively interact with inde- 
pendent environments, and since each qubit individually 
dephases in a Markovian way, the combined dynamics 
is also Markovian. For shorter qubit separations, the 
measure is non-zero and therefore the dynamics is not 
only indivisible (see Fig. [I] second row), but also non- 
Markovian in terms of back flow of information: more 
precisely, Nblp ^ whenever Nrhp ^ 0. We can there- 
fore identify the reservoir memory inducing the informa- 
tion backflow with the correlated noise and, in particular, 
with the cross-talk term. 

When the BEC interactions are stronger, i.e. for larger 
values of the scattering length a b , the local dephasing is 
always non-Markovian as the back action of the reser- 
voir on the system is non-negligible also when the qubits 
are far apart. In this case the BLP measure tends to a 
constant non-zero value at the point where environment 
mediated interactions vanish (D/L « 30). The value of 
non-Markovianity saturates and Mblp does not change 
for greater distances between the qubits. 



B. Non-Markovianity Dependence on Scattering 
Length 

In order to fully characterize the behavior of informa- 
tion flux, we now study how this quantity changes when, 
for fixed distances, we vary the scattering length of the 
ultracold reservoir gas. It is worth mentioning that while 
a change in the distance between the qubits can be seen as 
a change in the effective cutoff frequency of the reservoir 
spectrum, changing a b modifies the form of the spectrum 
itself and, in particular, its Ohmicity character [16] . It 
has been shown that for purely dephasing qubits the form 
of the spectrum, and more precisely its low-frequency 
component, has a strong effect on the presence of infor- 
mation backflow [5T]. Therefore we expect Nblp to be 
sensitive to changes in as, as can be already seen from 
Fig. 2. 

In Fig. [3] we plot the BLP measure as a function of the 
scattering length as for two values of qubits separations 
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FIG. 2. (Color online) Non-Markovianity Measure TVblp (the largest value in each column of the plot) as a function of distance 
D/L between qubits, for scattering lengths 0.02ag b < ttg"' (locally Markovian dynamics) and 0.50ag b > Bg"' (locally non- 
Mar kovian dynamics). The graph shows 20000 randomly drawn pairs of initial states. The pairs are classified into different 
categories and color coded accordingly as elaborated in the key on the right. 
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FIG. 3. (Color online) Non-Markovianity Measure TVblp (the largest value in each column of the plot) as a function of scattering 
length as between qubits for separation D = 200L (independent reservoir) and D = 4L (common reservoir). The graph shows 
20000 randomly drawn pairs of initial states. The pairs are classified into different categories and color coded accordingly as 
elaborated in the key on the right. 



D/L — 200 and D/L — 4, representing an independent 
(uncorrelated noise) and a common environment (corre- 
lated noise), respectively. For D/L = 200, if as > a% lt , 
we expect non-Markovian behaviour originating purely 
from the dynamics of the single qubit. There exists a 
cross-over between Markovian and non-Markovian pro- 
cesses, as found in Ref.[16], due to the single qubit dy- 
namical dependence on the scattering length. Only in 
presence of environment-mediated interactions between 
the two qubits (Fig. 3] , D/L = 4) the combined system 
can exhibit non-Markovian dynamics for < a c £ tl . In- 
deed, we see that even for very weakly interacting gases, 
we can have Nblp ^ 0. Hence, the presence of correlated 
noise makes the system non-Markovian for all values of 
the scattering length. 



Generally, for as > a^ lt , both single qubit 
non-Markovianity and environment-mediated non- 
Markovianity contribute to the total values of Nblp 
in a non-trivial way. While the single qubit local 
non-Markovianity, governed by the negativity of 71 (t), 
monotonically increases with a b after the threshold value 
corresponding to Markovian to non-Markovian crossover, 
non-Markovianity of the composite two qubits system in 
presence of correlated noise initially increases, reaches 
a maximum value and then decreases with as (Fig. [3j 
D/L = 200). This can be traced back to the behavior 
of the time dependent coefficients. In particular, for 
increasing scattering length, the environment-mediated 
interactions become weaker (72(0 is smaller) and hence 
the contribution of this term to Nblp decreases. 
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C. Optimizing pair 

In this section we examine the optimization procedure 
on which the BLP measure is based by studying numeri- 
cally the pair of states for which the backflow of informa- 
tion is maximized. This is arguably the most complicated 
system for which the optimizing pair has been studied in 
full detail. Our analysis shows that the optimizing pair 
is extremely sensitive to even small variations in the val- 
ues of the open system parameters. More specifically, 
the pair optimizing the BLP measure alternates between 
two different pairs of states when changing both the scat- 
tering length and the distance between the qubits. This 
underlines the practical difficulties in using the BLP mea- 
sure for increasing number of qubits. 

We have compelling numerical evidence that there are 
two candidates for the pair maximizing the measure. The 
first pair ($ + ,$~) corresponds to the sub-decoherent 
Bell states (decohering with rate r_(i)), while the sec- 
ond pair ( x & + , V E' _ ) comprises of the super-decoherent Bell 
states (decohering with rate r+(i)). The analytical ex- 
pression for the BLP measure is then 

A/blp = max{Af^ ,7V*}, (6) 

where A$ = ^ e - r -( b i) — e - r -( Q i) j s th e measure if the 
sub-decoherent Bell states form the maximizing pair and 
A* = Y^ri e -r +( hi ) — e _r+( - a ^ is the measure if the super- 
decoherent Bell states form the maximizing pair. Time 
intervals t £ indicate the periods of information 

backflow, manifested as dT±(t)/dt ~ 7i(£) ± 72 (0 > 0. 

We demonstrate the sensitivity of the maximizing pair 
on the model parameters in Fig. [4] by plotting the A/* 
and A/* as a function of the qubit separation. When 
the qubits share a common environment, corresponding 
to D/L < 30, the pair maximizing the measure depends 
strongly on the distance between the two qubits, alter- 
nating between the two sub-decoherent and the super- 
decoherent pair. This can be traced back to the dynam- 
ical behavior of 72 (t) and its absolute value in relation 
to the value of Ji(t). Due to the complex evolution of 
72 (t) the relative values of the two decoherence func- 
tions r± (/) oscillate with increasing distance, which is 
directly reflected in the relative values of A/* and A/*. 
For large qubit separations, the decoherence factors con- 
verge in the absence of environment-mediated interac- 
tions, r ± (t) -> 2F (i). Consequently for D/L > 30, the 
differences between A* and A/* diminish until in the case 
of effectively local environments both pairs maximize the 
measure. 

Expression ^ also demonstrates the qualitative 
agreement of the BLP and RHP measures of non- 
Markovianity. Recall that the latter is non-zero if either 
of the decay rates 7i(i)±7 2 (i) takes temporarily negative 
values. This immediately implies that either A$ or AAj> 
is non-zero, therefore proving that for the model consid- 
ered here non-divisibility of the two-qubit map exactly 
coincides with information backflow. 



D. Additivity issues for the non-Markovianity 
measure 

As a final step to complete the picture on two-qubit 
non-Markovianity we study the relationship between the 
degrees of non-Markovianity of the two-qubit channel, 
denoted A/2, and the local single qubit channel denoted 
by A/i . One of the ultimate goals of the theory of non- 
Markovian open systems is to contribute to better quan- 
tum technologies. A typical task in quantum information 
processing must eventually extend to larger multipartite 
systems consisting of n qubits and therefore it is of great 
importance to understand the scaling properties of non- 
Markovianity measures. When a given measure requires 
optimization over initial states, the task quickly becomes 
impractical for a larger number of qubits. Is it there- 
fore desirable to find links between the amount of non- 
Markovianity in multi-qubit channels and single qubit 
channels. 

When the two-qubit map factorizes in the limit of 
large qubit separation, = $a(£) ® &B(t), we dis- 

cover a simple connection between the single qubit non- 
Markovianity and the two-qubit non-Markovianity : 

A/2 = (e~ r ° (6 - ) + e _r °^)A/i. (7) 

Recall that To(t) is the single-qubit decoherence factor, 
and the corresponding decay rate has at most, a single 
period of negativity in the ultracold realization of the de- 
phasing model considered, i.e., 71 (t) < when t£ [a, b]. 
For as < a,% lt we have A/i = leading to A2 = 0: 
the two-qubit map is Markovian in the absence of single- 
qubit non-Markovianity and environment mediated inter- 
actions. When as > a C B %t and the single-qubit channel 
is non-Markovian, e _r °( fc ) + e~ r °( a ) < 2. As we can see 
then, even when the two-qubit channel factorizes into two 
identical single-qubit channels, the BLP measure is not 
additive, i.e, the information back-flow associated two 
independent qubits is not twice the back-flow of informa- 
tion associated with a single qubit. Instead, in this case, 
the measure is sub-additive, 2 A/i > A/2. 

Additivity properties of the map are even more com- 
plex in the common environment scenario, when the two- 
qubit map does not factorize. Interestingly, in this case 
the memory effects arising from correlated noise can lead 
to superadditivity: 2 N\ < N2 (see Fig. [5]). Moreover, 
the closer the qubits are to each other the larger is the 
difference |2 N\ — A/2], demonstrating the ability of the 
shared environment to enhance memory keeping prop- 
erties of the two-qubit map. It is an interesting open 
question if the sub-additivity of the BLP measure in 
the local environment scenario and super-additivity in 
the common environment case extends beyond the two- 
qubit case. An affirmative answer would imply that non- 
Markovianity scales more favorably in the common en- 
vironment case. Understanding the additivity properties 
of non-Markovianity measures is significant when using 
memory as a resource and our results highlight the im- 
portance of a suitable physical realization of a qubit reg- 
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FIG. 4. Maximising States: Non-Markovianity measure TVblp = max {A* , A/* } (the highest value at each instance) as a function 
of distance between qubits, D, for scattering lengths 0.02ag b < Og"' (locally Markovian dynamics) and 0.50ag b > ag' 1 (locally 
non-Markovian dynamics). Red line represents A/# and the blue line represents A/*. 




FIG. 5. Additivity of the non-Markovianity measure: A/* (red dotted line), A/* (blue dashed line) and 2A/i (black solid line) 
as a function of scattering length as, for qubit separations of D — 200L, D = 20L and D = 4L. Recall that A/2 = Ablp = 
max{A/$ , A* } . 



ister: the scaling of the non-Markovianity measure can 
crucially depend on whether the qubits share the same 
environment or not. 



V. CONCLUSIONS 

We have studied a system of two qubits dephasing in 
a structured environment. By explicitly taking into ac- 
count the spatial separation between the qubits we are 
able to interpolate between a common environment sce- 
nario, with the qubits close to each other, and an inde- 
pendent environment case, when the qubits are too far 
apart to influence each other. In the former case the dy- 
namics is characterized by correlated noise, arising from 
environment mediated correlations. Our focus was on 
the degree of non-Markovianity of the two-qubit map, 
quantified either in terms of non-divisibility of the map, 
or by the amount of information that the system can re- 
cover form the environment. These two characterizations 
of non-Markovianity, which generally do not agree, were 



shown to coincide qualitatively for the model considered 
here. 

Specifying our study to a physical realization of this 
model using optically trapped ultracold gases, we sin- 
gled out the microscopic origin of memory effects. In 
this case non-Markovianity emerges from the interplay 
of correlated noise and single qubit memory effects due 
to a non-trivial structure of the reservoir. When the lo- 
cal single-qubit dynamics is Markovian, two-qubit non- 
Markovianity can only emerge from correlated noise. 
More generally, the degree of non-Markovianity is am- 
plified when the qubits share an environment, highlight- 
ing the key role of correlated noise in non-Markovian 
systems. Moreover, the effect of the common environ- 
ment is also reflected in the global behaviour of the non- 
Markovianity measured based on information flux. Ac- 
tually, when the two qubits are very close to each other, 
the non-Markovianity measure is super-additive. Inter- 
estingly, when the qubits are in local environments the 
situation is reversed and the measure is sub-additive. Ad- 
ditivity properties affect the scaling of non-Markovianity 



measures for large qubit registers, potentially influencing 
the use on non-Markovianity as a resource for quantum 
technologies in a crucial way. We also discussed the pair 
of states optimizing the BLP measure and demonstrated 
the extreme sensitivity of the optimal pair on the system 
parameters. 
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explore the dynamics. The qubit takes physical form 
as an impurity atom occupying a double well potential 
with an optical supcrlatticc of wavelength A. The dis- 
tance between lattice sites, i.e the size of the qubit, is 
L = A/4. The minimum distance between qubits, where 
the two qubits occupy four consecutive lattice sites, is 
2D min = 8L. The superlattice is immersed in a Bosc- 
Einstein condensed (BEC) environment. We assume the 
condensate is trapped in a shallow potential and can be 
consider to be homogeneous. We also assume the gas 
to be weakly interacting, justifying the Bogoliubov ap- 
proach. The specific forms of 71 (i) and 72 (t) for this 
system, originally derived in Ref. |15j are 



Appendix A: Physical model 

In this work we consider a specific physical realization 
of the dephasing two-qubit model in order to explicitly 



7l(t) = ?kw r dkk 2 e -kv/2 MM^M (i_ M2kL)\ 



hir 2 J (e fe + 2g E n Q ) \ 2kL J 



M = 9se"o f°, H 2 - fc VV 2 sin(fr*)cos(|ft) / S in(2fc(£> + L)) sin(2fc(_D — L)) o sin(2fcJ) \ 
l2[ ) 2Tm 2 J (e k + 2g E n ) \ 2k(D + L) 2k(D — L) 2kD )' 



where g E — ^Ti 2 aE/rnE is the boson-boson coupling 
constant for a BEC environment with scattering length 
cle and mass m E , 9se — 2nh 2 asE/'msE is the coupling 
between the condensate (environment) and the impurity 
atom (system) with scattering length clse and reduced 
mass m S E = rn s m E /{ms + m E )- E k = ^2e k n g E + ef 



is the energy of the fc-th Bogoliubov mode, where no 
is the condensate density and e k = fi 2 k 2 / '(2mg). Fi- 
nally, a is the variance parameter of the lattice site. We 
consider specifically 23 Na impurity atoms immersed in a 
87 Rb condensate, with A = 600nm and no = 10 20 m~ 3 . 
The natural scattering length of the Rubidium atoms is 
a m = 99a where a is the Bohr radius, and we further 
assume cise — 55ao- 
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